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We investigate how a quenched random field influences the damage spreading transition in kinetic 
Ising models. To this end we generalize a recent master equation approach and derive an effective 
field theory for damage spreading in random field systems. This theory is applied to the Glauber 
Ising model with a bimodal random field distribution. We find that the random field influences 
the spreading transition by two different mechanisms with opposite effects. First, the random field 
favors the same particular direction of the spin variable at each site in both systems which reduces 
the damage. Second, the random field suppresses the magnetization which in turn tends to increase 
the damage. The competition between these two effects leads to a rich behavior. 



o 

CD 



CO 



I 

C 

O 

o 



> 

On 
On 
(N 

in 
o 

On 
c3 



-a 
c 

o 
o 



X 



05.40.+j, 64.60.Ht, 75.40.Gb 

The central question of damage spreading (DS) Jl)-^ is 
how a small perturbation in a cooperative system changes 
during the time evolution. This is analogous to the ques- 
tion to what extent the time evolution depends on the 
initial conditions, one of the main questions in non-linear 
dynamics that lead to the discovery of chaotic behavior 
||. In order to study DS the simultaneous time evolution 
of two replicas of a cooperative system is considered. The 
two replicas evolve stochastically under the same noise 
realization (i.e. the same random numbers are used in 
a Monte-Carlo procedure). The differences in the micro- 
scopic configurations of the two replicas are then used to 
characterize the dynamics and to distinguish regular and 
chaotic phases, depending on external parameters. 

Among the simplest cooperative systems are kinetic 
Ising models where DS has been investigated quite in- 
tensively within the last years using different dynamical 
algorithms such as Glauber |^,^-|| or heat-bath dynamics 
[^|,^|^Jl^] . In contrast to the equilibrium critical behavior 
the results of DS do depend on the particular choice of 
the dynamical algorithm although recently an attempt 
has been made to give a more objective definition of DS 
[fnj . In general, there are two different mechanisms by 
which damage can spread in a kinetic Ising model. First, 
the damage can spread within a single ergodic compo- 
nent (i.e. a pure state or free energy valley) of the system. 
This is the case for Glauber or Metropolis dynamics. Sec- 
ond, the damage can spread when the system selects one 
of the free energy valleys at random after a quench from 
high temperatures to below the equilibrium critical tem- 
perature. This is the only mechanism to produce DS in 
an Ising model with heat-bath dynamics. This algorithm 
is thus well suited for exploring the structure of the free 
energy landscape. 

In the literature the name DS has been applied not 
only to the studies discussed above but also to a differ- 
ent though related type of investigations in which the two 
systems are not identical. Instead, one or several spins in 
one of the copies are permanently fixed in one direction. 
Thus the equilibrium properties of the two replicas devi- 
ate from each other and their microscopic differences can 



be related to equilibrium correlation functions |12| , [13[ . 
Note that in these works the use of identical noise (i.e. 
random numbers) for the two systems is not essential but 
only a method to reduce the statistical error. 

Whereas DS in clean Ising models is comparatively well 
understood less is known about disordered models. The 
influence of random fields has been investigated in a two- 
dimensional Ising-like model with Metropolis dynamics 
giving a reduction of the damage at high temperatures 
but an increase at low temperatures p] . By using the 
heat-bath algorithm DS has been used to study the phase 
space structure of Ising spin glasses [jl5|— |lT|| and the cor- 
responding critical behavior at the DS transition [jl8j . 

In this Letter we consider the original DS problem, viz. 
the time evolution of two identical systems and study 
the influence of a quenched random field on DS in ki- 
netic Ising models. To this end we generalize the master 
equation approach |?],|| to random field systems. The 
resulting effective field theory of DS is then applied to 
the Glauber Ising model with a bimodal random field 
distribution. We study the dependence of the spreading 
transition on temperature and field and determine the 
phase diagram. 

We consider two identical Ising models with N sites 
described by the Hamiltonians and H^ 1 given by 
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where S} is an Ising variable with the values ±1 and 
n = f,2 distinguishes the two replicas. Jy is the (non- 
random) exchange interaction between the spins. The 
random field values (pi are chosen independently from a 
distribution p[ip). The dynamics of the systems are given 



by stochastic maps S, 
Glauber algorithm 



(n) 



(i + f) = F[{S 7 (rl) (*)}],e.g., the 
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where the transition probability v(x) is given by the usual 
Glauber expression 



v(x) 



Here h["\t) — J2j JijS)' LI {t) + ipi is the local magnetic 
field at site i and (discretized) time t in the system n. 
£i(i) £ [0,1) is a random number which is identical for 
both systems, and T denotes the temperature. 

Within the master equation approach |?],[| the simul- 
taneous time evolution of the two replicas is described by 
the probability distribution 
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where (•} denotes the average over the noise realizations. 
The variable with the values H — h, H — , — h, or de- 
scribes the states of the spin pair (SV,SP). Thedis- 
tribution P fulfills the master equation 



dt 



P(vi, ...,v N ,t) = 



N 



-/] } J P(yi, ■ ■ ■ ,Vi, ■ ■ ■ ,VN,t)w(Vi —> Hi) 
N 

+ ^2 ^2 P( v iT--^iT--,VN,t)w{H t ->u i ). (5) 

i — 1 [JLi^Vi 

The transition probabilities w(fii — > Vi) have to be calcu- 
lated from the properties of the stochastic map F which 
defines the dynamics. 

As in the clean case we derive an effective field theory 
by assuming that fluctuations at different sites are statis- 
tically independent which amounts to approximating the 
distribution P{v\, ■ ..,ujf,t) by a product of single-site 
distributions P Ui (i,t). However, in a disordered system 
different sites are not equivalent and thus their P Vi (i,t) 
are not identical. This is the main difference to the clean 
case where the single-site distributions are all iden- 
tical. In the following we further assume that P Vi {i,t) 
is determined by the local value tpi of the random field 
only, P Vi [i, t) = P Vi (ipi, t). In general, this is an approxi- 
mation since sites with identical random field values may 
well have different environments which should influence 
the distribution. In the mean-field limit of infinite dimen- 
sions oder infinite-range interactions, however, the above 
assumption becomes exact. 

Inserting the decomposition 



(6) 



i=l 



into the master equation (||) gives a system of coupled 
equations of motion for the single-site distributions 



+ P ll (<p)W(n^v,<p)], (7) 

where W(fi — ► ip) is the transition probability w aver- 
aged over the states z/, of all sites. We now define the 
local damage d(ip) = P_| (ip) + P ^(<p). The total dam- 
age is obtained as the corresponding disorder average 



= / d(pp(<p)d(<p). (8) 




The equation of motion of the local damage can be easil 
determined from (ffl). Using some symmetry relations 
between the transition probabilities W , it reads 



dt 
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So far the considerations have been rather general, be- 
ing valid for any (single-site) dynamic rule and any dis- 
tribution of the random field. We now apply the formal- 
ism to the Glauber Ising model in the mean-field limit 
Jij = Jo/N (for all i and j). The random field distri- 
bution remains unspecified so far. In order to determine 
the spreading point for infinitesimal initial damage it is 
sufficient to solve ([}]) in linear order in d(ip). After cal- 
culating the transition probabilities W in analogy to the 
clean case and some further algebra we obtain 

j t d(<p) = -\m(<p)\ d(ip) + ^[1 - m 2 (^)] D. (10) 

If we concentrate on DS processes starting in equilibrium 
conditions the local magnetization 



m(ip) = tanh[(Jom + f)/T] 



(11) 



and the average magnetization m are time- independent. 
Eq. (1lC|) is very similar to the corresponding eq. (36) of 
Rcf. jpffor DS in a homogeneous field. The main differ- 
ence is that for random-field systems we have to distin- 
guish between the local damage d(<p) which determines 
the healing probability [first term in dlQ)] and the average 
damage D which determines the damaging probability of 
a site [second term in (|l0|)]. In a homogeneous system 
local and average damage are identical. Consequently, 
replacing d(ip) by D and m(ip) by m in (|^|) exactly gives 
the corresponding equation for the homogeneous system. 

To proceed further we have to specify the random field 
distribution p(ip). As an example we will discuss the 
bimodal distribution 



p(<f) = ~ <Po) + S{<P + <Po)} fao>0). (12) 

The thermodynamics of the mean-field Ising model 
with a bimodal random field has been investigated in 
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FIG. 1. Thermodynamic phase diagram of the mean-field 
Ising model with bimodal random field. TCP denotes the 
tricritical point. 



FIG. 2. Lyapunov exponents of the mean-field Glauber 
Ising model with bimodal random field. The peak in the curve 
for ifio = 0.4 corresponds to the T where m_ vanishes. 



detail almost 20 years ago |1J 
takes the form 



The equation of state 
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The resulting phase diagram is summarized in Fig. [l]. 
There is a tricritical point at Ttcp = 2 Jo/3 and ^tcp ~ 
0.439Jo- For T > Ttcp the ferromagnetic phase transi- 
tion is of second order, for T < Ttcp it is of first order. 

We now turn to our results on DS in this model. Using 
the notations d± — d(±tpo)i d = and m± — 

m(±(/3 ) the equation of motion (M) can be written as 



dt 



A d. 



(14) 



The dynamical matrix is given by 



A = 



(l-m 2 _)^ -\m-\ + (1 -m 2 _) 



2T 



(15) 



The question whether the damage spreads or heals can 
be answered by means of the eigenvalues of A. If both 
eigenvalues are negative the damage heals, if at least one 
of them is positive the damage spreads. 

In the paramagnetic phase we have |m+| = |m_| = 
tanh(<^ /T). The eigenvalues of A are given by Ai = 
— m+ + (1 — m 2 ^) Jo/T and A2 = —rn+. The correspond- 
ing eigenmodes are the average damage D and the dam- 
age difference d + — cL, respectively. Consequently, in 
the paramagnetic phase the Lyapunov exponent which is 
given by the largest eigenvalue of A reads 

A = - tanh^o/T) + [1 - tanh 2 (</> / T )] J /T. (16) 

Its dependence on temperature and random field strength 
is visualized in Fig. ||. In the paramagnetic phase 



the Lyapunov exponent decreases with increasing ran- 
dom field ipo and therefore the spreading temperature T s 
which is defined by A = increases. This can be easily 
understood from the fact that a random field favors a 
particular spin direction at each site. Since this direction 
is the same for the two replicas the corresponding spins 
in the two replicas tend to be parallel which reduces the 
damage. For random field strength ipo > Jo the Lya- 
punov exponent remains negative for all temperatures 
and thus the damage never spreads. Asymptotically for 
<Po — * Jo we obtain 



Tl = - J? 



1 
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(17) 



We note that the critical tpo which completely suppresses 
DS has the same value as the corresponding critical ho- 
mogeneous field H although the functional dependence 
of T s on the field is different. 

In the ferromagnetic phase \m+ \ and |m_ | are different. 
In this case the eigenvalues of A are given by 

Ai,2 = ^ - |m_| + ^(1 -m+ + 1 -m 2 _) 

± \\(\m+\ \m-\) 2 + ^(1 - m 2 + + 1 - m 2 _) 2 



- |m_|)(m+ - ml) 



1/2 



(18) 



In order to calculate the Lyapunov exponent we first de- 
termine the average magnetization m as a function of <po 
and T from the equation of state ( |l3| ) . We then calculate 
to + and 77i— and insert them into (|l^). The resulting 
Lyapunov exponents are presented in Fig. ^. In contrast 
to the paramagnetic phase the spreading temperature de- 
creases with increasing random field strength. At a first 
glance this seems to contradict the argument given above, 
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FIG. 3. Damage spreading phase diagram of the mean-field 
Glauber Ising model with bimodal random field 



viz. that a random field favors a particular spin direction 
and thus reduces the damage. However, the random field 
also influences DS via a reduction of the magnetization 
since the Lyapunov exponent (jl^) is determined by the 
local magnetizations. In the ferromagnetic phase this 
effect is stronger than that of the preferred orientation 
discussed above and thus T, is reduced. 

we have determined the 



By means of (16) and 
spreading temperature as a function of the random field 
strength. The resulting phase diagram of DS in the 
mean-field Glauber Ising model with bimodal random 
field is shown in Fig. ||. The minimum spreading tem- 
perature Tg^jnin « 0.438 is obtained when the spreading 
transition coincides with the ferromagnetic phase transi- 
tion which occurs for <p w 0.480 (see Fig. §. 

To summarize, we investigated the influence of a 
quenched random field on DS in kinetic Ising models. We 
generalized the master equation approach to random 
field systems and derived an effective field theory for DS. 
As an example we studied the mean-field Glauber Ising 
model with bimodal random field. We found that the 
random field supports the spreading of damage in the 
ferromagnetic phase but hinders it in the paramagnetic 
phase. For strong enough field the damage never spreads. 

In the concluding paragraph we discuss other random 
field distributions and compare our results to the numer- 
ical simulation |lj] . The influence of the particular form 
of the random field distribution on DS can be discussed 
qualitatively by means of . This equation shows that 
the healing probability is proportional to the local mag- 
netization. This means that the damage on sites with 
local magnetization zero cannot heal. Consequently, DS 
will be qualitatively different in systems with a continu- 
ous random field distribution since even for very strong 
random fields there will be sites with vanishing local mag- 
netization. Thus damage will spread on a subset of sites 
with low enough random field. However, with T — > 
the measure of this subset goes to zero. A detailed in- 



vestigation of this case will be published elsewhere p0[ . 
These results also help to understand the numerical sim- 
ulation which was carried out for a box distribution. 
It shows a decrease of the spreading temperature with 
increasing random field although the stationary value of 
D is reduced at high temperatures. This is consistent 
with a reduction of T s due to a suppression of the local 
magnetization and spreading on a subset of sites at low 
temperatures. However, a direct comparison with the 
mean-field theory is not possible since in the simulations 
a two-dimensional system was used which - due to fluc- 
tuations - does not have an ordered phase for any finite 
random field. Finally, we discuss possible extensions of 
this work. Besides a systematic investigation of different 
random field distributions the damage equation of motion 
should be solved beyond first order in the damage. This 
will permit the determination of the stationary damage 
values and the investigation of the critical behavior at 
the spreading transition. Some studies along these lines 
are in progress Q . 
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